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ABSTRACT 

We present the morphological properties of 161 star clusters in M33 using the Ad- 



CZ5 

C$ ' vanced Camera For Surveys Wide Field Channel onboard the Hubble Space Tele- 

scope using observations with the F606W and F814W filters. We obtain, for the first 
^ ■ time, ellipticities, position angles, and surface brightness profiles for a significant 

QO ■ number of clusters. On average, M33 clusters are more flattened than those of the 

| Milky Way and M31, and more similar to clusters in the Small Magellanic Cloud. 

I/"") , The ellipticities do not show any correlation with age or mass, suggesting that rota- 

O 

^vq . tion is not the main cause of elongation in the M33 clusters. The position angles of 

the clusters show a bimodality with a strong peak perpendicular to the position angle 
of the galaxy major axis. These results support the notion that tidal forces are the 
$_j \ reason for the cluster flattening. We fit King and EFF models to the surface bright- 



ness profiles and derive structural parameters including core radii, concentration, 
half-light radii and central surface brightness for both filters. The surface brightness 
profiles of a significant number of clusters show irregularities such as bumps and 
dips. Young clusters (Log age < 8) are notably better fitted by models with no ra- 
dial truncation (EFF models), while older clusters show no significant differences 
between King or EFF fits. M33 star clusters seem to have smaller sizes, smaller 
concentrations, and smaller central surface brightness as compared to clusters in the 
MW, M31, LMC and SMC. Analysis of the structural parameters presents a age- 
radius relation also detected in other star cluster systems. The overall analysis shows 
differences in the structural evolution between the M33 cluster system and cluster 
systems in nearby galaxies. These differences could have been caused by the strong 
differences in these various environments. 
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1 INTRODUCTION 

The spatial structures and internal stellar kinematics of star clusters (SCs) bear the imprint of their 
initial conditions and their dynamical evolution. Determination of structural parameters of SCs 
yields important information to understand the processes involved in their formation and evolu- 
tion. Due to their proximity, galaxies in the Local Group provide us with ideal targets for detailed 
studies of SC structural parameters. While the SC systems of the Milky Way (MW) and M31 have 
received close attention, the third s piral galaxy in the Local Group, M33, has been less studied 



As po inted out by previous authors (|Sarajedini et al 



1998 



Chandar et al 



2006; 



San Roman et al 



20101) . clusters in M33 populate almost all evolutionary phases, including ages and metallicities 
not observed in the MW or M31. Therefore, M33 star clusters can provide us with essential infor- 
mation to understand the formation and evolution of the ir early and late stages. At a distance of 



870 Kpc (distance modulus = 24.69; 



Galleti et al. 



2004), M33 is the only nearby late-type spiral 



galaxy (Scd) and provides us with an ideal opportunity to study the star cluster system of such a 
galaxy. 

The structure of SCs can be investigated through their surface brightness profiles. These pro- 
files are often characterized by two families of models. The first set, kno wn as King models, are 



based on tidally limited profiles of isothermal spheres (King 



1962, 



1966J). More extensive King- 



based models exist that take into consideration non-isotropic systems, rotating systems or even 



multi-mass models (e.g. 



Michie 



1963 



Wilson 



1975 



Da Costa & Freeman 



1976). The second fam- 



ily of models, known as EFF (for Elson, Freeman, and Fall) models, were empiricall y derived to 



reproduce the surf ace brightness profiles of young clusters in the Magellanic Clouds (|Elson et al. 



19871 : 



Elson 



19911) . These profiles are described as power-laws that do not include the pronounced 



tidal truncation of the King models. 

It is well known that King models provide an excellent des cription of the luminosity pro- 



files a nd internal kinema tics of most MW globular clusters (e.g. iDjorgovski 



1995 



1995 : 



Trager et al. 



McLau ghlinl |2000|) . High spatial resolution data are required to study spatial structures in 



extragalactic clusters. Hubble Space Telescope (HST) imaging has been successsfully used to fit a 
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variet y of models to the surface brightness profiles of clusters in M31 (IBarmbv et al 



2002 



2007 



2009|V the Magellanic Clouds (MC), and the F ornax dwarf spheroidal galaxy (|Mackey & Gilmore 



2003a b 



ck lMcLaughlin & van der Marelll2005h . These studies suggest that globular clusters (GCs) 



describe a fundamental plane analogous to the fundamental plane of elliptical galaxies. Morpho- 
logical studies of younger clusters present more difficulties based on their, generally, lower con- 
centration and mass * it is not well established if open clusters (P Cs) in the MW also lie on the same 



fundamental plane (IBonatto & Bica 



2005 



Barmbv et al. 



20091) . Therefore, it is essential to deter- 



mine the structural parameters of star clusters in other galaxies to constrain SC initial conditions 
and evolution. 

In the las t few years, significant advances have been made to u nderstand the star cluster system 



of M33 (e.g. 



San Roman et al 



this system are very limited. 



2009 



2010; 



Larsen et al. 



Cockcroft et al 



201 lb . However, structural studies of 



(|2002|) present structural parameters derived from King 



model fits for four GCs in M33. Dynamically and structurally, t hey appear virtually id entical to the 



Chandar et al 



(|1999|) present core 



MW and M3 1 GCs, and fit very well into the fundamental plane, 
radii of 60 star clusters using linear correlations with the measured full width at half maximum 
(FWHM) of each cluster. For these reasons, we have undertaken a comprehensive structural study 
of th e largest sample of M33 star clusters to date. Throughout this work, a distance to M33 of 870 



Kpc (IGalleti et al.ll2004) has been adopted. At that distance, 1" corresponds to 4.22 pc. 



This paper presents the results of this study and is organized as follows. Section 2 describes the 
observations and data reduction. Section 3 presents the analysis of ellipticities and position angles 
while Section 4 discusses the surface brightness profiles and the profile fittings. The analysis of the 
structural parameters and comparison with other galaxies are shown in Section 5. Finally, Section 
6 presents a summary of the study. 



2 OBSERVATIONS AND DATA REDUCTION 

he observations and d ata reduction for the present study are described in a companion paper 



(San Roman et al 



2009, hereafter Paper I), and for a detailed description we refer the reader to 
that paper. For convenience, we provide an abbreviated summary below. 

The observations were obtained with the Advanced Camera for Surveys Wide Field Chan- 
nel (ACS/WFC) on board the HST Twelve HST/ ACS fields from the GO- 10190 program (RL: 
D. Garnett) have been analyzed. Three filters were used for the primary observations (F475W, 
F606W, F814W) and two for the parallel images (F606W, F814W). The fields are located along 
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the southwest direction covering a radial extension between 0.9 - 6.6 Kpc from the center of the 
galaxy. Figure 1 in Paper I shows the locations of these fields while Table 1 of that paper presents 
a summary of the observations. 

All of the images were calibrated through the standard pipeline-process. The standard-calibrated 
"FLT" images were corrected for the geometric distortion and then used to obtain the point- 
spread function (PSF) . A detailed description of the PSF construction method is presented in 



Sarajedini et al.l (|2006|) . This PSF was used to convolve the structural models before the fitting 



process. 

The pipeline-processed drizzled ("DRZ") images were used to construct the surface brightness 
photometry. We derived positional offsets between these frames using the imshift and imcombine 
tasks in IRAF to allow us to produce one master image per filter per field. To make the calculation 
of the photometric errors more straightforward, each DRZ image was multiplied by the exposure 
time before performing the photometry, and then the background sky value was added back to 
each image. 



3 ELLIPTICITIES 



Crowding and incompleteness in the inner region of the clusters make star counts impractical. 
For this reason, we have derived the structural parameters of the clusters from luminosity profiles 
rather than number density profiles. A precise determination of the cluster centers is crucial to 
obtain this morphological information of the objects. Inaccurate centers would produce artificial 
distortions in the radial profiles. To obtain cluster centers and shapes we used the IRAF task EL- 
LIPSE to fit elliptical isophotes to our images in both filters. ELLIPSE is an algorithm designed to 
fit galaxy profiles that decrease monotonically, so we first filtered the images with an 1 1 x 1 1 pixel 2 
median filter to create the required smooth profiles. Spatial smoothing of the images is a common 
technique when dealing with surface brightness rather than number density of stars, however the 
required magnitude of the smoothing could degrade the accuracy of the derived centers. At the 
distance of M33, our objects are partially unresolved so our images do not need a high degree 
of filtering to obtain an appropriate s mooth profiles. This technique has been analyzed and sue 



cessfu lly used in previous studies (e.g. 
2008) . 



Larsen et al 



2002; 



Barmbv et al 



2002; 



McLaughlin et al. 



Isophotes were fitted between 0".l and the largest measurable semi-major axis in both bands. 
The overall center position of each cluster has been determined as the average of the ELLIPSE 
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Figure 1. Ellipticities and position angles as a function of semi-major axis for a representative subsample of the total cluster sample. Solid lines 
correspond to values from F606W images while dotted lines correspond to F814W images. Position angles are counted north through east. The 
identification name of each cluster correspond with the notation in Paper I. 



output between 0".5 and 1".5 along the semi-major axis where the measurements are more reliable. 
The final averaged position of each cluster differs by less than 0". 1 between the F606W and F8 14W 
filters. Because of this small offset in the cluster position between the different filters, we have 
chosen the cluster coordinates determined from the F606W filter for the following analysis. 

A second pass of ELLIPSE was run to determine ellipticities and position angles (PA), this 
time fixing each cluster center to the previously calculated values. The ellipticity is a complex 
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quantity that can vary at different radii from the cluster center. To analyze any potential variation, 
elliptical isophotes were fitted between 0".l and the largest measurable semi-major axis distance 
in both filters. Figure \T\ shows the ellipticity and PA profiles as a function of semi-major axis for 
a representative sample of clusters. The arbitrary fluctuations in the ellipticities and PAs at small 
radii (< 0".5) are likely to be produced by internal errors in the ELLIPSE algorithm. In the far outer 
parts of the clusters, the ellipticity and PA are poorly constrained due to the low signal-to-noise 
ratio. Given these facts, the final ellipticity and PA for each cluster was calculated as the average 
of the values between 0".5 and 1".5 where the quantities are more stable. In most of the cases, the 
ellipticity is well traced in the different filters. In 12 cases, the faintness of the objects prohibited us 
from an accurate determination of the ellipticity in both the F606W and F814W filters. For small 
ellipticities, the PA is not well determined and varies significantly between filters. We chose not 
to use PA values which had errors exceeding 20%, and hence we have excluded those PA values 
from the final catalog. Table[T|shows the ellipticities and PAs for the star cluster sample identified 
in Paper I in the F606W and F814W filters. Errors correspond with the standard deviation of the 
mean. 



Some studies have suggested t hat cluster flattening is 



effects rather than galactic tides ([White & Shawl 



1987 



generally attributable to cluster rotation 



Davoust & Prugniel 



99pL I 



dynamical evolution of these objects, they lose mass and angular momentum. As a consequence, 
clusters rotate slower and become rounder as they evolve. To investigate the elongation origin of 
M33 SCs, Figure [2] shows different cluster properties as a function of ellipticity for F606W. The 
cluster properties were determined and analyzed in Paper I. Ellipticities in M33 clusters do not 
seem to be driven by any of the cluster properties plotted in Figure |2] 



The distributions of cluster flattening in nearby galaxies are known to be very different (e.g. 
van den Berghll2008r) . Ellipticities in Galactic GCs have been found to correlat e inversely with age , 



suggesting structural changes over time or changes in the initial conditions (jFrenk & Fall 



1982) 



The correlat ion of the ellipt i city with the age and mass in LMC and SMC clusters are found 



to be weak ( Goodwin 



1997 



Hill & Zaritskv 



2006b . Our M33 sample does not show any clear 



correlation with age or mass. A lumi nosity - ellipticity relation has also been suggested in some 
nearby galaxies (Ivan den Berghll2008|) . While the faintest clusters associated with massive galaxies 
(MW, M31 and NGC 5128) are flatter than the most luminous star clusters, no evidence for a 
luminosity - ellipticity correlation exists in the LMC or SMC. No correlation between magnitude 
and ellipticity is found in our sample either. Elsonl ( 1991 ) suggest that the high ellipticity observed 
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Figure 2. Cluster ellipticities versus mass (M), age, galactocentric distance (R gc ), V magnitude and V-I color. The cluster properties were adopted 
from Paper I. 



in many of the young LMC clusters could be due to the presence of subclumps. However, such 
substructures will be erased as clusters evolve. 

The mean ellipticity of star clusters in the MW/M3 1 is smaller than that of the LMC, which 
is in turn smaller than the mean ellipticity of star cl usters in the SMC. It has been suggested tha t 



this is correlated with the masses of these galaxies (|Geisler & Hodge 



1980; 



Staneva et al 



Goodwin 



1996) 



19971) suggests the tidal field strength of the parent galaxy as the dominant factor driving 



the differences between the LMC and the MW star clusters. Tidal forces would destroy the initial 
anisotropic velocity dispersions of a cluster, modifying its initial shape. This effect would be much 
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Figure 3. Mean ellipticity of star clusters as a function of their host galaxy luminosity. 



less significant in the less massive galaxies so the tidal field in the LMC would not be able to 
modify significantly the shapes of its clusters maintaining their initial triaxiality. This scenario 
would also explain the even higher ellipticities observed in the SMC. Figure |3] shows the mean 
ellipticity of different star cluster systems as a functio n of the lumi nosity of their host galaxy . 



The mea n ellipticities w e re obt ained from 



Harris 



for M31, 



van den Bergh 



1996 ) for the M W, 



Hill & Zaritsky 



Barm by et al. 



(2002, 



(|2006h for the SMC. The mean 



(I2008h for the LMC and 
ellipticity for clusters in M33 corresponds to the mean value from Table 1 from F606W band. 
The error of the data correspond s with the standard deviation. The luminosities for each galaxy 
are adopted from Ivan den Berghl (|1999|) . With a mean ellipticity of e=0.2, M33 clusters are, on 
average, more flattened than those of the MW or M31 and more similar to clusters in the MC in 
agreement with this scenario. 

On the other hand, tidal fields are also be able to stretch clusters and make them more elon- 
gated. If tidal forces have a dominant effect on the elongation of M33 clusters, one would expect 
a preferred orientation of the SC position angles towards the location of the force. Figure |4] sum- 
marizes the elongation and orientation of the clusters with respect to the galaxy. There is no clear 
trend in the orientation vectors toward the galactic center, however the orientation of the clusters 

© 201 1 RAS, MNRAS 000.Q]-?? 



Newly Identified Star Clusters In M33. Ill 9 



i — i — i — i — i — i — i — i — i — | — i — i — i — i — i — i — i — i — i — | — i — i — i — i — i — i — i — i — i — | — i — i — i — i — i — i — i — i — r 




C_ i i i i i i i i i I i i i i i i i i i I i i i i i i i i i I i i i i i i i i i 

23.6 23.5 23.4 23.3 23.2 

Right Ascension (degrees) 

Figure 4. Cluster elongations and orientations represented in the RA-Dec coordinate plane. North is up and east is to the left. The vector sizes are 
correlated with the ellipti city obtained in this st udy. The cross corresponds with the center of M33, with the arrow marking the position angle (23°) 
of the galaxy major axis jRegan&Vogelll994l) . 



tend to have a preferred direction. In particular, the position angle distribution in Figure |5] shows a 
bimodality with a strong peak at - 55". This northwest direction points toward M31, which is lo- 
cated at a position angle of about - 40° with respect to the center of M33 . Recent evidence suggests 



that a close encounter between M33 anc 



(Putmanetal 



2009; 



McConnachie et al 



M31 could have affected the properties of the M33 disk 



20091 ) . Due to the small galactocentric distance of these 



clusters (R gc < 7 Kpc), the idea of M31 as the source of the elongation seems unlikely. Analysis of 
the PAs show no correlation with age, mass, distance or magnitude of the clusters. No correlation 
has been found with the spiral arm pattern or with circular orbits around the galaxy. The elongation 
of the clusters are not aligned with the x-y pixel coordinate system of each set of images, ruling 
out a possible artifact related to cluster orientation. The significance of the preferred orientation of 
the M33 clusters remains open and further investigation of this phenomenon will require precise 
PAs of a bigger sample of clusters at different locations in the galaxy. 
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Figure 5. Distribution of the position angles (PA) of our sample. 



4 SURFACE BRIGHTNESS PROFILES 

As mentioned before, we have derived surface brightness profiles rather than surface number den- 
sity profiles. The ELLIPSE task was again run to obtain F606W and F814W surface brightness 
profiles. This time the fixed, zero-ellipticity isophote mode was required be cause we chose to fit 



Barmbv et al. 



(2007) to transform the raw 



circular models. We have followed the prescription in 
output values from ELLIPSE to calibrated surface brightnesses on the Vegamag system. We have 
first converted image counts (counts pixeL 1 ) into luminosity density in L Q pc~ 2 . We have worked 
in terms of linear intensity and finally transform the results to surface brightness magnitude units 
(///mag arcsec" 2 ). The profiles extended out to R=14" (63pc) in both filters. Figure |6] shows the 
cluster surface brightness profiles where the error bars are estimations obtained by ELLIPSE from 
the isophotal intensities. 



4.1 Irregular Profiles 

Visual inspection of the surface brightness profiles shows that a significant number of cluster 
profiles do not decrease smoothly as expected, exhibiting irregularities such as bumps, dips or 
sharp edges. These irregularities appear in both filters. Figure|6]shows examples of more traditional 
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profiles versus Figure [7J that shows examples of cluster profiles with different irregularities. For a 
detailed inspection of the irregualr profiles, we have overp lotted 3 concentric apertures, 0.1", 0.5" 
and 1 .0", on the F606W images in Figure |7J More than half of our sample exhibits some level 
of anomaly. Some of these features could be the result of statistical fluctuations or due to a few 
luminous stars. In fact the oversampling of our data could magnify some features in the profiles, 
conferring an occasional peak due to a very bright star, the degree of a pronounced bump. However, 
sub-clumps or a deficit of stars are also prominent in a significant number of clusters (e.g. clusters 
12, 24, 71), and appear to reflect real structures within the clusters. 



Si milar anomalies have been identified among LMC a nd SMC star cluster profiles (|Mackey & Gilmore 



2003b; 



Hill & Zaritskv 



2006 



Werchan & Zaritskv 



20111) . They detect, in several surface bright- 



ness profiles, systematic deviations from the analytical models that they attribute to a lack of 
central concentration. They refer to these clusters as "ring" clusters and suggest that these anoma- 
lies could be the result of dynamical evolution. Some of the clusters could be in the process of 
dissolution or have had insufficient time to dynamically rela x. Vari a tions i n the mass to light ra- 



Elson 



(1199 1|) reports bumps, sharp 



tio could also cause departures from the analytical models, 
"shoulders", and central dips in the profiles of LMC clusters with an age range similar than our 
sample and interpret these anomalies as consequences of the initial conditions or signatures of 
merging sub-condensations. As the clusters evolve, the substructure will be erased, but we do not 
see any correlation between the presence of these anomalies and any other parameters. In par- 
ticular, objects with irregularities are observed along the entire range of ages with no preference 
at any particular stage of evolution. Whether any physical significance can be attached to these 
irregularities is uncertain. 



4.2 Profile Fitting: King vs EFF 



Several different mo dels h a ve bee n used to descri be star cluste r density profiles. The mor e com 



mon choices include 



hereafter Wilson), Elson et al. 



Kind( 



19621. he reafter King). lKingl ( 



1987 



1966J, hereaf ter King66), 



Sersic 



Wilson! (11975 



1968, hereafter Sersic). 



hereafter EFF) and 

The empirical King model has been frequently used to describe the number surface density 
profiles of old globular clusters. King66 and Wilson models are single-mass, isotropic, modi- 
fied isothermal sphere, defined by a phase-space distribution function. Wilson model differs from 
King66 in an additional term that produces clusters spatially more extended than King66 model. At 
small and intermediate radii, the two model profiles are quite similar. King66 and Wilson models 
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Figure 6. Surface brightness profiles and model fits to a representative sample of clusters. Left panels correspond with F606W profiles and right 
panels with F814W. The solid line in each panel represents the best-fit King model while the dotted line represents the best-fit EFF model. 



have n o an alytic expressions for den sity or surface brightness. See McLaughlin & van der Marel 



(2005) and 



McLaughlin et alJ (|2008r) for a detailed comparison of both models. The other two 
models, EFF and Sersic, correspond with "power-law" models. They are defined to observation- 
ally describe more extended surface-density profiles. Like the Wilson model, Sersic model was 
originally designed for galaxy application. 

The family o f King model (King66 and K i ng) ar e the most commonly used in star clusters 



studies, however 



McLaughlin & van der Marell (|2005|) find that Wilson models fit clusters of any 
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Figure 6. Continued 



age, in any galaxy, as well as or better than King66 models. McLaughlin et alJ (|2008r) conclude 
that analyses of structural correlations will not be significantly affected by using King66 models 
when Wilson profiles are better fits or vice versa. Based on this discussion we decide to fit a tidally 
truncated model and a more spatially extended model. For simplicity we chose the empirical King 
and EFF models to describe th e surface bri ghtness profiles of our sample. 



The empirical King model (IKing 



1962 ) is described as: 
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Figure 7. Same as Figure[6]but for a group of clusters showing irregularities in their profiles. Each row corresponds with a different cluster with 
left panel showing the F606W profile and right panel the F814W profile. The bottom row shows a 3"x 3" F606W image of the 3 clusters. Each 
image is shown with the same gray-scale with north up and east to the left. For visual guiding, the 3 concentric apertures correspond with radii of 
0.1", 0.5" and 1".0. 
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P(R)=Po-{ r rl (1) 

{[l+(R/R c ) 2 Y 2 [l+(R,/Rc) 2 ]' 2 ) 

where p is the central surface number density of the cluster, R c is the core radius, and R, is the 
tidal radius. It is useful to define the concentration parameter as c = log(R t /R c ). 

However, young LMC cluste rs do not appear to be tidally truncated, and seem to be more suit- 
ably fitted by power-law profiles. Elson et al.l (|1987|) argue in favor of an unbound "halo" caused by 
the cluster expansion due to mass loss or violent relaxation, and proposed an empirical description 
known as the EFF model: 

fi(R)= M0 -(i+^y (2) 

where fio is the central surface brightness, a is the scale radius, and y is the power-law slope 
at large radii. The discrepancies between the two profiles are significant at large radii, being es- 
sentially indistinguishable in the inner regions. The derived half-light radius, Rf„ can be estimated 
based on the transformation equation in Larsen 2006. For the King model: 

R h = 0.541 -R c • (iW) - 486 (3) 
and for the EFF profile: 

R h = a • V(0.5)^ - 1 (4) 

Considering the presence of young clusters in our sample we decided to fit the cluster profiles 
with both types of models via weighted x 2 minimization. An additional constant term, <f> bk , was 
added to the models to compensate for the background contamination. Although the clusters are 
quite well resolved, the effect of the point-spread function (PSF) in the core regions cannot be 
ignored. To impr ove the accuracy of the fitted parameters, the models were convolved with an 



appropriate PSF (|San Roman et al.l 120091 for details) before the minimization process. From the 
total sample of 161 clusters in Paper I, 33 objects were rejected due to low signal-to-noise ratios. 
Figure |6] shows the cluster surface brightness profiles with the best-fitting models where the error 
bars are estimations obtained by ELLIPSE from the isophotal intensities. Clusters with atypical 
profiles, Figure [7] deviate from the analytical models in specific regions. These clusters have, 
in general, slightly larger x 2 values, relative to both models, than clusters with no anomalies. In 
fact the fit to some of these clusters are almost certainly not physically descriptive but even in 
these cases, the fits are sufficiently good to characterize the sizes of the clusters and their cores. 
The parameters from the best-fit King model are given in Table El with Table [3] presenting the 
parameters from the best-fit EFF model. We report the calibrated and extinction-corrected central 
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Figure 8. Quality of the fit of EFF model relative to King model versus age and half-light radius for the two different filters. 



surface brightness in the native band-pass of the data. Reddening values have been adopted from 
Paper I. When indi vidual reddening is no t available, we adopt a line-of-sight reddening value 
of E(V - I) = 0.06 (lSarajedini et al.ll2000|) . The listed errors in the tables, are given by the x 2 
minimization fitting process. Alternatively, the difference between the F606W and F814W fits can 
be used as a more realistic uncertainty. 

Nine clusters, with available structural parameters, are common to 



Chandar et al 



(1999). Di 



rect comparison between core radii gives us offsets of: A/? C (CBF - Us) = - 0.81+ 0.3 pc. The pro 



Chandar et al. 



files of those objects, except in two cases, have prominent irregularities in their profiles. 
(|l999 ) calculate core radii using a linear correlation with FWHM rather than direct model fits. The 
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Figure 9. Comparison of half-light radius observed in both F606W and F814W filters for the different fitted models. 



discrepancy between the two studies can be caused by the different adopted methods in addition 
to the peculia r profiles. 



Following McLaughlin & van der Marell (|2005|) . we define a statistic parameter that compares 
the^ 2 of the best fit of both models for any object: 



2 2 

Xeff ~Xk 



nig 



Xeff + Xk 



(5) 



nig 



If the parameter is zero, the two models fit the same cluster equally well. Positive values in- 
dicate a better fit of the King model, and negative values indicate EFF models are more suitable. 
Figure [8] shows this A parameter as a function of age and half-light radius. Estimated ages for our 
M33 sample were adopted from Paper I, where they were obtained via color-magnitude diagram 
(CMD) fitting. Young clusters (Log age < 8) are notably better fit by models with no radial trunca- 
tion. In older clusters, small values of the statistic parameter, -0.2 < A < 0.2, show no significant 
differences between the quality of the fittings. EFF models are also favored for smaller sizes as 
shown in the bottom panel of Figure [8] 

Working with F606W and F814W provides us with two independent sets of measurements 
for each cluster. Systematic errors and color dependencies in the fits can be evaluated through 
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direct comparison of both measurements. Figure |9] compares the parameters derived from model 
fits in F606W and F814W. The results show, in general, good agreement between bands, however 
specific clusters with significant disagreement between filters stand out. These differences in the 
half-light radii are more prominent for the EFF derived parameters. In the majority of the cases 
correspond with clusters with strong irregular profiles in both bands or anomalies in one filter 
and not in the other. In particular the derived structural parameters of clusters 3, 28, 33, 34, 70, 
71, 95, 109, 116, 124, 159 are significantly in disagreement for the two different bands. We only 
consider the F606W model fits from now on because this filter is more sensitive to the underlying 
stellar population, and less sensitive to foreground contaminating stars. In addition, fits to clusters 
in other galaxies are typically performed in the V-band allowing us direct comparison without 
being concerned about possible color gradients. Since King models are primarily used to describe 
clusters in other galaxies, we adopt the derived King structural parameters as the basis for the 
subsequent analysis. 



5 STRUCTURAL PARAMETERS 



5.1 Comparison with Other Galaxies 

The dynamical evolution of star clusters directly affects their structural properties. Therefore, com- 
paring structural parameters in different galaxies is of special importance to understand star cluster 
evolution in a variety of environments. We have compared the newly derived M33 structural pa- 



rameters with those for clusters in M3 1 ( Barmbv et al 



20071) . the LMC and SMC (IMcLaughlin & van der Marel 



20051) . M31 sample has been analyzed in an homogeneous way in the same filters F606W and 
F814W. Comparisons with LMC and SMC samples correspond with direct comparison with V- 
band photometric quantities. In addition, this study is based in number density converted to net 
mag arcsec 2 . Both studies have fitted a number of different models including King66. 

While King model is described by the core radius, R c , King66 is characterized by the scaled 
radius R . It is important to recognize that R and R c are not strictly equivalent. For the King model, 
the core radius corresponds with the radius at which the density has fallen to just half the central 
density, lKm g (Rc) = 0.5 • Iq, while for King66 this relation is only approximate, lKing66(Ro) ~ 0.5 • I . 



The theoretical models agree closely with t he empirica 



values and fit the observations equally well (|King 



family of curves for a reasonable range of 



1966). For this reason R is often called the core 



radius. The connection between the theoretical and the observational parameters is straightforward 
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to calculate and both studies of comparison provide, in addition to the basic King66 fit parameters, 
the derived structural parameters. 

Figure [10] shows correlations among diffe r ent str uctural parameters. We have included the four 



old, halo M33 GCs analyzed in 



Larsen et al. 



(120021) . They appear to have similar structural prop- 
erties than our sample except for the higher brightness. All the galaxies present similar trends 
although several differences stand out. M33 seems to possess clusters with small sizes, not show- 
ing objects with radii between R c = 3 - 10 pc observed in other galaxies. With an average of c= 1.1 2, 
the concentration in M33 clusters is smaller than the mean concentrations of clusters in the MW 
(c=l .41), M31 (c=l .46) and LMC (c=l .46). The top-right panel in Figure H0]shows that the central 
surface brightness of M33 star clusters is fainter than the one in M3 1 and LMC star clusters, cover- 
ing a smaller range and presenting a slightly offset position in the diagram. Direct comparison with 
V-band photometric quantities, in the case of LMC and SMC samples, could motivate the offset. 
In general more massive/luminous clusters tend to have brighter /u . The brighte st clusters in M33 



are m ore than 1 mag fainter than the brightest clusters in MW, M31 and LMC (ISan Roman et al. 



2010|) so one would expect M33 star clusters to have fainter central surface brightness. 



Structural properties for the M33 sample are shown as a function of estimated age in Figure 
ITTl When possibl e, we have compared our sample with the above studies. Ages for M31 clusters 
are adopted from Caldwell et al.l (|2009l) for the common objects, while ages for LMC and SMC 
are provided in their previously mentioned studies. 

Comparison between these galaxies sugg ests similar trends . The age-radius relation among 



M33 clusters has been previously sugge st ed (IChandar et al 



galaxies (e.g. 



Mackev & Gilmore 



2003aUbl : 



Glatt et al. 



1999 ) and also observed in different 



2009|), suggesting that the age - radius trend 



Elsonetal 



( 1989 ) iden- 



has its origin in internal cluster processes, rather than external elements, 
tified, among LMC clusters, a clear trend that t he spread in core radius is a increasing function 
of cluster age. Thi s trend was confirmed later by DVlackey & Gilmorel (|2003a|) using HST images. 



Elson et al 



(|1989t) suggests mass loss from stellar evolution as the cause of cluster expansion. All 



clusters would have formed wi th small core radii an d as the most massive stars evolve they would 



expand. In a more recent study, 



Mackev et al 



(|2008|) use N-body simulations to specifically model 
the evolution of massive star clusters observed in LMC and SMC. They identify two physical pro- 
cesses which result in substantial and prolonged core expansion. The first process, mass-loss due 
to rapid stellar evolution in a primordially mass-segregated cluster, only occurs during the 100 - 
200 Myrs after the formation of the cluster. The second process, heating due to a retained popula- 
tion of stellar mass black holes, does not begin until several hundred Myr after the formation but 
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Figu re 10. Structural parameters of the present study as compared with star clusters in M3 1 jBarmbv et aD l2007n . LMC and SMC 
iMcLaughlin & vim der Marel200^) . R c and R/, are plotted in logarithmic scale for better comparison. 



result in core expansion for many Gyr. A combination effect of these two process result in a wide 
variety of evolutionary paths that cover the observed distribution of massive star clusters in LMC 
and SM C. Several other mechanisms have been tested as possible explanations. IWilkinson et i 
( 2OO3I) investigate the effect of a time- varying tidal field and variations in primordial binary frac- 
tions in N-body simulations. Assuming a point-mass galaxy the clusters show similar core radius 
evolution on both circular and elliptical orbits and therefore they conclude that the tidal field of the 
LMC does not influence the evolution of the clusters significantly. The presence of large numbers 
of primordial binaries in a cluster leads to core radius expansion; however the magnitude of the 
effect is insufficient to explain the observations in LMC clusters. 

The trend of older clusters having a larger range in core radii than the younger population is 
clearly visible in the top-right panel of Figure [TT] The youngest clusters in our sample possess 
compact cores of R c ~ - 1 pc. At intermediate ages M33 clusters expand a range of < R c < 4, 
slightly larger increase than among the LMC clusters in the same age period. A faster mass-loss 
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Figu re 11. Structural parameters of the present study versus age as compared with M31 (Bar mbvetafl |2007|) . LMC and SMC 
iMcLaughlin & vim der MarefcobH) 



rate could produce a rapid evolution of M 33 clusters. Variations in the initial mass function (IMF) 
slope could be responsible for this effect. iGieles & Baumgardtl (|2008|) . in an extensive survey of 
N-body simulations, show that the strength of the tidal field is the dominant factor determining the 
escape rate of stars. This result suggests that external elements can also produce an enhancement 
in the mass-loss rate and accelerate the evolution of the cluster. Strong environmental effect could 
also affect the evolution of M33 star clusters. 



5.2 Galactocentric Distribution 

Star clusters lie on the tidal field of their host galaxy and therefore are subject to tidal shocks. 
It is then reasonable to expect that their position from the center of the host galaxy can influ- 
ence their dynamical evolution. Clusters closer to the center of the galaxy will be more dynami- 
cally evolved due to the effects of the tidal field. Structural parameters of MW globular clusters 
have been shown to be independent of galactocentric distance except for the half-light radius (e.g. 
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van den Bergh et al 



1991 



galaxies (|Barmby et al. 



2002, 



Diqrgovski & Meylanlll994|) . Similar results have been found in nearby 



20071) . These studies show that star clusters at large galactocentric 



distances, on average, have larger diameters than those closer to the galactic center. While large 
clusters with small galactocentric distances could have been disrupted, there is no satisfactory ex- 
planation for the deficiency of small clusters at larger distances. For this reason, this correlation 
has been suggested to be an initial physical condi tion of cluster formation and a general property 



of star cluster systems (|van den Bergh et al. 



1991). 



Only projected distances are available for our sample. No correlation has been found between 



concentration and galactocentric distance, R gc . Centr al surface brig^ 



ltness does not correlate either 



Barmbv et al. 



(|2007l) . and use the renormal- 



with R gc . To explore the R^ - R HC relation, we follow 
ized galactocentric distance as a better indicator of galactic potential than galactocentric distance 
itself: 

R gc /(SKpc) 



D* 

gc ~ VJillOkms" 1 ) 



(6) 



where V c corresponds to the galaxy's c ircular velocity. We have u sed the followi ng values of 
(inkms" 1 ): 220 for the MW 100forM33 (Corbelli & SalucciLoOO ). 23 for M3 1 ( Carignan et al 



20061) . 65 for the LMC (Ivan der Marel et al.ll2002l) and 60 for the SMC (IStanimirovic et al. 



2004). 



Figure d2] shows R/, as a function of renormalized galactocentric distance, R* c , for the different 
galaxies. The filled circles in the top left panel represe nt our sample of M3 3 clusters, and the star 



Larsen et al 



(I2002|) . The solid line corre- 



symbols correspond with the four globular clusters in 
sponds with the least-squares fit of the combined sampl e. The solid lines for the rest of the galax- 



ies have been taken from Table 10 (Barmbv et al. 



2007|) as the suggested star cluster fundamental 



plane fits. No clear correlation has been found in the M3 3 sample. Only 5 confirmed clusters have 



been detected in the o utskirts of M33 (R ;r > 10 Kpc) ( 



Cockcroft et al. 



2011). 



Cockcroft et al. 



Stonkute et al 



2008 



Huxor et al 



2009; 



(|2011|) fit King models to these clusters and derived struc- 
tural parameters. When these 5 clusters are included in the analysis, the data show a similar trend 
than in other galaxies (Top left panel of Figure [12]). The five outer clusters are represented as open 
symbols, and the new least-squares fit corresponds with the dashed line. While the trend in other 
galaxies is clearly driven by small inner clusters and also large clusters at larger radii, the corre- 
lation in M33 is just forced by very large clusters in the outer part of the galaxy. This could be a 
selection effect since at large distances, small clusters will be difficult to detect in ground-based 
images. We note that the star cluster systems in this analysis correspond with different age ranges, 
however as Figure QT] shows there is no correlation between R h with age. Therefore any differ- 
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Figure 12. Half-light radius as a function of renormalized galactocentric distance for clusters in our sample, M31, MW and MC. Solid lines 
represe nt the least-squares fits. In the top left panel, the solid line represents the least-squares fit including our sampl e (filled circles) and th e four 
GCs in Larsen et al. (2002)(star symbols). Open circles in the same panel correspond with outskirt clusters in M33 l lCockcroft et al.ll201 ll) . The 
dashed line corresponds with the least-square fit including the outskirt cluster in M33. 



ences bet ween galaxies in Fig ure [T2l are not due to secondary correlation with this parameter. In 



addition, 



Barmbv et al. 



(2002, 



2007 ) find a mild correlation between R h and [Fe/H] in the MW, 



M31 and the MC. iBrodie & Straderl (|2006|) discuss several possible explanations and argue that 
this correlation is a consequence of two effects, the - R gc correlation and the predisposition of 
metal rich clusters to occupied smaller galactocentric radius. 

High-resolution imaging has shown that the overall properties of GCs in different galaxies are 
remarkably similar. In fact, several studies suggest that GC structural and dynamical properties 



can be summarizec 



(|Djorgovski 



1995 



as describing a fundamen t al plane analogous to the one for elliptical galaxies 



Dubath & Grillmair 



the four M33 clusters studied in 



1997 



Larsen et al 



McLaughlin 



20001) . Dynamically and structurally, 



(|2002|) . appear virtually identical to MW and M31 



globular clusters. Measurements of velocity dispersions for a large sample of M33 clusters would 
permit precise M/L ratios to determine with confidence if M33 clusters fall in the suggested fun- 
damental plane for MW and M31 clusters. 
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6 SUMMARY 

We present ellipticities, position angles and structural parameters of M33 star clusters using HST/ACS 
images in F606W and F814W bands. This study presents for the first time detailed morphologi- 
cal properties of a significant number of M33 star clusters. Surface brightness profiles have been 
fitted with King models and EFF models to determine structural properties such as as core and 
tidal radii, concentration and central surface brightness. On average, the clusters in M33 are more 
elliptical than those of the MW or M31 and more similar to clusters in the SMC. No correlation 
was found between the ellipticities and any of the other properties studied herein, including age 
or luminosity. The position angles of the clusters show a bimodality with a strong preferred ori- 
entation perpendicular to the position angle of M33 major axis. A significant group of clusters 
shows irregularities such as bumps, dips or sharp edges in their surface brightness profiles. Young 
clusters are notably better fitted by models with no radial truncation (EFF models), while older 
clusters show no significant differences between King or EFF fits. M33 star clusters seem to have 
smaller sizes, smaller concentrations, and smaller central surface brightness as compared to clus- 
ters in the MW, M31, LMC and SMC. The sample presents an age-radius relation also detected 
in other galaxies. Structural parameters show no correlation with galactocentric distance although 
a trend with the half-light radius appears when combining our sample with M33 outer clusters. 
The overall analysis shows differences in the structural evolution between the M33 cluster system 
and cluster systems in nearby galaxies. These differences could have been caused by the strong 
differences in these various environments. 
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Table 1: Ellipticities and Position Angle (PA). 



Id a 


Ellipticity 


Ellipticity 


PA 


PA 




(F606W) 


(F814W) 


(F606W) 


(F814W) 


1 


n o o i n n i 

0.23 ± 0.01 


n on i n no 

0.29 ± 0.02 


63.4 ± 2.6 


64.2 ± 2.2 


2 


A O C i f\ f\ A 

0.25 ± 0.04 








3 
4 


n 1 a ■ a ao 

0.14 ± 0.03 


n 1 o in n a 
[J. lis ± 0.04 






5 
6 


A 1 *7 ■ n 

0.17 ± 0.02 


a 1 ~7 n no 

0.17 ± 0.03 






7 


f\ OO i n AO 

0.23 ± 0.02 


n o c n n i 

0.25 ± 0.01 


64.5 ± 2.0 


62.7 ±1.5 


X 


n f\c i a a 1 

0.05 ± 0.01 


n 1 /i in no 

0.14 ± 0.02 






9 
10 


A 1 n , A AO 

0.17 ± 0.03 








1 1 


0.14 ± (J. 02 


a on i a no 

0.20 ± 0.02 






12 


n on i n f\A 

0.30 ± 0.04 


n o o n no 

0.33 ± 0.03 






13 


n 1 /" in no 

0.16 ± 0.03 








14 


A h-\/^ , a AO 

0.22 ± 0.03 


n oo n no 

0.22 ± 0.02 


63.1 ±4.1 


62.4 ±5.4 


15 


n ao i n no 
0.43 ± 0.03 


n o o i n no 

0.38 ± 0.02 


A A n i A O 

-40. y ±4.3 


-52.7 ± 5.7 


16 


n o i n no 

0.21 ± 0.02 


A OO . A AO 

0.22 ± 0.03 






17 


n 1 c in no 
0.15 ± 0.02 


a 1 o in no 

O.lo ± 0.03 






1 8 


n on i n n 1 

0.29 ± 0.04 








19 


n i o in no 

0.13 ± 0.02 


A 1 A i A A1 

0.10 ± 0.01 


-46.3 ±3.9 


-45.6 ± 3.8 


20 


n 1 n i f\ r\A 

0. 19 ± 0.04 


n n a no 

0. 17 ± 0.03 






2 1 


n o j^t i n no 

0.36 ± 0.02 








22 


n r\£. i n m 

0.06 ± 0.01 


A AA i A AO 

0.09 ± 0.02 






23 


a no . n no 

0.28 ± 0.03 








24 


n o a in no 

0.24 ± 0.02 


a 1 c i n no 

0.15 ± 0.02 






25 


n i n in no 

0. 19 ± 0.02 


n ot n no 

0.25 ± 0.03 


41.7 ± 4.2 


37.8 ± 5.6 


26 


a , r\ no 

0.26 ± 0.03 


n on i n no 

0.29 ± 0.02 


-62.0 ± 4.3 


-58.2 ± 4.4 


27 


n i o in no 
0.12 ± 0.03 


n 1 o in no 
0.13 ± 0.02 






28 


r\ i a in nn 
0.14 ± 0.00 




- /0.0 ± 0.0 




29 


n o c i n nc 

0.25 ± 0.05 


n o o i n r\A 

0.33 ± 0.04 






30 


T\ + 01 


1% + 01 


-69.4 ± 3.2 


-66.4 ± 3.6 


31 


0.24 ± 0.03 


0.26 ± 0.02 






32 


0.24 ± 0.02 


0.26 ± 0.02 






33 


0.09 ± 0.01 




-30.0 ±5.1 




34 










35 


0.27 ± 0.03 


0.37 ± 0.02 






36 


0.14 ±0.00 


0.22 ± 0.00 


-70.0 ± 0.0 


-75.1 ±0.0 


37 


0.23 ± 0.02 


0.18 ± 0.03 


67.5 ± 1.2 


-28.0 ± 25.2 


38 


0.12 ±0.00 


0.17 ± 0.03 


-58.3 ± 0.0 


-23.6 ± 17.4 


39 


0.45 ± 0.02 


0.40 ± 0.01 


38.6 ±3.2 


30.9 ± 0.4 


40 


0.21 ± 0.02 


0.19 ±0.03 






Note. 


- Table[T]is published in its entirety in the electronic edition of the Journal. 



Identification number in Paper I 
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Table 2: Structural Parameters from the King Profile Fit. 



to 





Filter 






R 






R, 




c 




R;, 


R;, 


<Pbk 






(mag arsec~~) 


(arcsec) 


(pc) 


(arcsec) 


(PC) 








(arcsec) 


(pc) 


(L pc- 2 ) 


1 


F606W 


19.17 


±0.02 


0.22 ± 0.00 


0.92 ± 0.02 


2.47 


+ 


0.09 


10.43 ± 0.36 


1.06 


+ 


0.02 


0.39 ± 


0.01 


1.63 ±0.03 


34.96 ± 0.25 




F814W 


18 55 


+ 02 


0.21 ± 0.00 


0.90 ± 0.02 


2 41 




OS 
u.uo 


10.16 ±0.33 


1.05 


+ 


0.02 


U.JO :£ 


01 


1.60 ±0.03 


39.82 ± 0.28 


2 


F606W 
F814W 


1 8 97 
18 89 


+ 02 
+ 01 


0.10 ±0.00 
0.14 ±0.00 


0.41 ± 0.01 
0.57 ± 0.01 






















37.47 ± 0.28 
42.16 ±0.26 


3 


F606W 
F814W 


21 19 

20 78 


+ 06 
+ 10 


0.79 ± 0.04 
0.32 ± 0.02 


3.35 ±0.16 
1.34 ±0.08 


3 84 




19 


16.19 ± 0.79 


0.68 


+ 


0.03 


93 + 


03 


3.94 ±0.13 


36.90 ± 0.32 
42.50 ± 0.38 


4 


F606W 


18.08 


+ 0.02 


0.09 ± 0.00 


0.39 ± 0.01 


2 83 




22 


11.92 ±0.92 


1.48 


± 


0.04 


27 + 


01 


1.13 ±0.04 


36.71 ±0.20 




F814W 


18 15 


+ 01 


0.10 ±0.00 


0.44 ± 0.01 


3 00 




37 


12.65 ± 1.55 


1.46 


+ 


0.05 


29 + 


02 


1.23 ±0.07 


43.17 ±0.26 


5 


F606W 


20 43 


+ 03 


0.19 ±0.00 


0.79 ± 0.02 






















42.60 ± 0.36 




F814W 


20 43 


+ 08 


0.22 ± 0.01 


0.94 ± 0.05 


4 00 




1.65 


16.86 ±6.95 


1.25 


+ 


0.18 


50 + 


10 


2.09 ± 0.42 


50.11 ±0.39 


6 


F606W 


19 21 


+ 02 


0.22 ± 0.00 


0.94 ± 0.02 


4.60 




0.67 


19.41 ± 2.84 


1.32 


+ 


0.06 


53 + 


04 


2.23 ±0.16 


40.15 ±0.29 




F814W 


1 8 83 


+ 03 


0.22 ± 0.01 


0.91 ±0.03 


7 37 




2 64 


31.09 ± 11.14 


1.53 


± 


0.16 


65 + 


11 


2.76 ± 0.48 


45.54 ± 0.39 


7 


F606W 
F814W 


21.66 
21 72 


+ 07 
+ 06 


0.38 ± 0.02 
0.38 ± 0.01 


1.60 ±0.08 
1.58 ±0.06 






















25.56 ± 0.22 
29.47 ±0.17 


8 


F606W 


19.66 


+ 03 


0.29 ± 0.01 


1.21 ±0.03 


2 22 




0.08 


9.38 ± 0.35 


0.89 


± 


0.02 


0.42 + 


0.01 


1.79 ±0.04 


39.05 ± 0.38 




F814W 


19 17 


+ 03 


0.27 ± 0.01 


1.13 ±0.03 


1 95 




0.06 


8.24 ± 0.27 


0.86 


± 


0.02 


39 + 


01 


1.63 ±0.03 


43.73 ± 0.43 


10 


F606W 


20 03 


+ 04 


0.49 ± 0.02 


2.05 ± 0.07 


6 19 




79 


26.10 ±3.33 


1.11 


± 


0.06 


92 + 


06 


3.86 ± 0.25 


41.73 ±0.25 




F814W 


20 22 


+ 05 


0.65 ± 0.02 


2.73 ±0.10 


5 11 




51 


21.57 ± 2.15 


0.90 


± 


0.05 


97 + 


05 


4.08 ± 0.21 


46.59 ± 0.34 


11 


F606W 


19.86 


+ 0.03 


0.61 ±0.01 


2.58 ± 0.06 


4 50 




22 


18.96 ±0.91 


0.87 


± 


0.02 


88 + 


02 


3.72 ±0.10 


44.13 ±0.39 




F814W 


19 17 


+ 03 


0.59 ± 0.02 


2.49 ± 0.07 


4 24 




20 


17.88 ±0.84 


0.86 


± 


0.02 


84 + 


02 


3.55 ± 0.09 


49.36 ± 0.57 


12 


F606W 


20 14 


+ 02 


0.22 ± 0.00 


0.91 ±0.02 


2 00 




0.06 


8.42 ± 0.27 


0.97 


± 


0.02 


35 + 


01 


1.47 ±0.03 


25.01 ±0.15 




F814W 


19.70 


±0.02 


0.22 ± 0.00 


0.95 ± 0.02 


1.83 


± 


0.05 


7.74 ± 0.22 


0.91 


± 


0.02 


0.34 ± 


0.01 


1.44 ±0.03 


28.48 ±0.17 


13 


F606W 


18.77 


± 0.03 


0.21 ± 0.01 


0.90 ± 0.03 


2.28 


± 


0.08 


9.60 ± 0.32 


1.03 


± 


0.02 


0.37 ± 


0.01 


1.56 ±0.03 


40.74 ± 0.38 




F814W 


18.38 


±0.02 


0.22 ± 0.00 


0.93 ± 0.02 


2.02 


± 


0.07 


8.50 ± 0.29 


0.96 


± 


0.02 


0.35 ± 


0.01 


1.49 ±0.03 


46.60 ± 0.50 


14 


F606W 


18.80 


±0.02 


0.44 ± 0.01 


1.84 ±0.04 


5.03 


± 


0.22 


21.20 ±0.94 


1.06 


± 


0.02 


0.78 ± 


0.02 


3.30 ± 0.08 


24.96 ± 0.29 




F814W 


19.17 


± 0.03 


0.48 ± 0.01 


2.01 ± 0.05 


5.85 


± 


0.43 


24.67 ± 1.81 


1.09 


± 


0.03 


0.88 ± 


0.03 


3.72 ±0.14 


183.60 ± 0.29 


17 


F606W 


20.71 


±0.03 


0.54 ± 0.02 


2.29 ± 0.08 


3.99 


± 


0.22 


16.85 ± 0.94 


0.87 


± 


0.03 


0.78 ± 


0.03 


3.30 ±0.11 


23.81 ±0.22 




F814W 






0.70 ± 0.03 


2.95 ±0.12 


3.31 


± 


0.08 


13.98 ± 0.34 


0.68 


± 


0.02 


0.81 ± 


0.02 


3.44 ± 0.08 


181.31 ±0.30 


18 


F606W 


18.23 


±0.03 


0.15 ± 0.01 


0.64 ± 0.03 


2.46 


± 


0.10 


10.39 ± 0.42 


1.21 


± 


0.03 


0.32 ± 


0.01 


1.35 ±0.04 


27.36 ± 0.41 




F814W 


18.68 


± 0.03 


0.18 ± 0.01 


0.76 ± 0.03 


2.08 


± 


0.06 


8.77 ± 0.27 


1.06 


± 


0.02 


0.32 ± 


0.01 


1.36 ±0.04 


185.03 ±0.40 


19 


F606W 


18.51 


±0.01 


0.19 ± 0.00 


0.79 ± 0.02 


4.66 


± 


0.45 


19.63 ± 1.91 


1.39 


± 


0.04 


0.49 ± 


0.02 


2.06 ±0.10 


45.80 ± 0.28 




F814W 


18.03 


±0.02 


0.20 ± 0.00 


0.86 ± 0.02 


4.03 


± 


0.33 


16.98 ± 1.38 


1.30 


± 


0.04 


0.47 ± 


0.02 


2.00 ± 0.08 


51.76 ±0.36 


22 


F606W 


18.72 


±0.01 


0.51 ± 0.00 


2.16 ±0.02 


5.46 


± 


0.17 


23.05 ± 0.70 


1.03 


± 


0.01 


0.89 ± 


0.01 


3.74 ± 0.06 


45.29 ± 0.28 




F814W 


17.97 


± 0.01 


0.54 ± 0.01 


2.27 ± 0.03 


5.31 


± 


0.15 


22.39 ± 0.63 


0.99 


± 


0.01 


0.89 ± 


0.01 


3.78 ± 0.06 


52.70 ± 0.39 



r 



Note. - Table[2]is published in its entirety in the electronic edition of the Journal. 
a Identification number in Paper I 



Table 3: Structural Parameters from the EFF Profile Fit. 





Filter )jlq 


a c 




y 


R/, 


R/i 


<f>bk 




(mag arsec -2 ) 


(arcsec) 


(pc) 




(arcsec) 


(pc) 


(L pc- 2 ) 



1 


F606W 


19.43 


+ 


0.02 


0.39 


+ 


0.03 


1.66 


+ 


0.1 1 


4.56 ± 0.38 


0.33 


+ 


0.04 


1.00 


+ 


0.00 


36.76 ± 0.22 




F814W 


18.80 


+ 


0.02 


0.35 


+ 


0.02 


1.48 


+ 


0.09 


4.16 ± 0.30 


0.33 


+ 


0.04 


1.00 


+ 


0.00 


41.90 ± 0.25 


2 


F606W 


18.93 


+ 


0.02 


0.08 


± 


0.00 


0.32 


+ 


0.02 


1.65 ± 0.06 


0.08 


+ 


0.00 


0.00 


+ 


0.00 


36.74 ± 0.29 




F814W 


18.78 


± 


0.02 


0.07 


+ 


0.00 


0.29 


± 


0.02 


1.29 ± 0.04 


0.06 


± 


0.00 


0.00 


+ 


0.00 


39.99 ± 0.40 


3 


F606W 


22.05 


± 


0.05 


1.19 


+ 


0.23 


5.00 


± 


0.98 


5.13 ± 1.47 


0.89 


+ 


0.29 


3.00 


+ 


1.00 


37.89 ± 0.26 




F814W 


20.85 


± 


0.04 


0.28 


+ 


0.03 


1.18 


± 


0.14 


1.81 ± 0.21 


0.28 


+ 


0.03 


1.00 


+ 


0.00 


41.91 ± 0.42 


4 


F606W 


18.20 


± 


0.02 


0.13 


± 


0.01 


0.54 


± 


0.03 


2.76 ± 0.10 


0.29 


+ 


0.05 


1.00 


+ 


0.00 


37.22 ± 0.19 




F814W 


18.27 


± 


0.02 


0.14 


+ 


0.01 


0.61 


± 


0.03 


2.79 ± 0.12 


0.31 


+ 


0.06 


1.00 


+ 


0.00 


43.67 ± 0.26 


5 


F606W 


20.38 


± 


0.03 


0.12 


± 


0.01 


0.51 


± 


0.04 


1.31 ± 0.08 


0.11 


± 


0.01 


0.00 


+ 


0.00 


41.20 ± 0.40 




F814W 


20.59 


± 


0.04 


0.21 


± 


0.03 


0.88 


± 


0.10 


2.06 ± 0.25 














49.67 ± 0.50 


6 


F606W 


19.33 


± 


0.02 


0.25 


+ 


0.01 


1.08 


± 


0.05 


2.50 ± 0.13 


0.98 


± 


0.41 


4.00 


+ 


1.00 


40.16 ± 0.32 




F814W 


18.90 


± 


0.03 


0.23 


± 


0.01 


0.97 


± 


0.06 


2.23 ± 0.12 














45.23 ± 0.41 


7 


F606W 


21.63 


± 


0.05 


0.32 


+ 


0.03 


1.37 


± 


0.11 


1.77 ± 0.12 


0.32 


± 


0.03 


1.00 


+ 


0.00 


25.28 ± 0.19 




F814W 


21.62 


± 


0.05 


0.27 


+ 


0.02 


1.15 


± 


0.10 


1.47 ±0.10 


0.26 


± 


0.02 


1.00 


+ 


0.00 


28.93 ± 0.22 


8 


F606W 


20.02 


± 


0.02 


0.56 


+ 


0.06 


2.38 


± 


0.27 


6.04 ± 0.98 


0.36 


+ 


0.06 


1.00 


+ 


0.00 


41.32 ± 0.26 




F814W 


19.56 


± 


0.03 


0.55 


± 


0.07 


2.33 


± 


0.31 


6.55 ± 1.28 


0.33 


± 


0.06 


1.00 


+ 


0.00 


46.55 ± 0.30 


10 


F606W 


20.25 


± 


0.02 


0.69 


± 


0.05 


2.89 


± 


0.23 


3.36 ± 0.30 


0.91 


± 


0.19 


3.00 


+ 


0.00 


42.04 ± 0.27 




F814W 


20.66 


± 


0.01 


1.15 


± 


0.17 


4.84 


± 


0.73 


5.41 ± 1.23 


0.81 


± 


0.20 


3.00 


+ 


0.00 


47.47 ± 0.35 


1 1 


F606W 


20.26 


± 


0.01 


1.34 


± 


0.13 


5.64 


± 


0.53 


7.24 ± 1.05 


0.74 


+ 


0.09 


3.00 


+ 


0.00 


46.08 ± 0.31 




F814W 


19.58 


± 


0.01 


1.39 


± 


0.15 


5.86 


± 


0.63 


8.19 ± 1.39 


0.70 


+ 


0.09 


2.00 


+ 


0.00 


52.48 ± 0.44 


12 


F606W 


20.45 


± 


0.01 


0.49 


± 


0.06 


2.08 


± 


0.25 


6.54 ± 1.08 


0.29 


± 


0.05 


1.00 


+ 


0.00 


26.12 ±0.10 




F814W 


20.07 


± 


0.01 


0.63 


± 


0.11 


2.67 


± 


0.45 


9.34 ± 2.43 


0.29 


+ 


0.06 


1.00 


+ 


0.00 


29.79 ±0.12 


13 


F606W 


19.02 


± 


0.03 


0.33 


± 


0.02 


1.41 


± 


0.10 


3.99 ±0.32 


0.34 


± 


0.04 


1.00 


+ 


0.00 


43.26 ± 0.30 




F814W 


18.66 


± 


0.02 


0.37 


± 


0.03 


1.57 


± 


0.13 


4.60 ± 0.50 


0.31 


± 


0.05 


1.00 


+ 


0.00 


49.83 ± 0.40 


14 


F606W 


19.01 


± 


0.02 


0.79 


± 


0.06 


3.32 


± 


0.25 


4.42 ± 0.39 


0.69 


± 


0.09 


2.00 


+ 


0.00 


26.03 ± 0.30 




F814W 


19.47 


± 


0.01 


0.81 


± 


0.06 


3.40 


± 


0.25 


4.13 ±0.36 


0.77 


+ 


0.10 


3.00 


+ 


0.00 


184.38 ±0.31 


17 


F606W 
F814W 


21.12 


± 


0.02 


0.96 


± 


0.10 


4.06 


± 


0.44 


5.37 ±0.85 


0.69 


+ 


0.12 


2.00 


+ 


0.00 


24.60 ± 0.20 


18 


F606W 


18.43 


± 


0.03 


0.50 


± 


0.05 


2.12 


± 


0.19 


6.77 ± 0.75 


0.29 


± 


0.03 


1.00 


+ 


0.00 


29.73 ± 0.26 




F814W 


19.01 


± 


0.03 


0.71 


± 


0.09 


3.00 


± 


0.36 


10.53 ± 1.92 


0.30 


± 


0.04 


1.00 


+ 


0.00 


187.80 ±0.27 


19 


F606W 


18.63 


± 


0.01 


0.25 


± 


0.01 


1.04 


+ 


0.04 


2.71 ±0.09 


0.60 


± 


0.11 


2.00 


+ 


0.00 


46.14 ±0.30 




F814W 


18.18 


± 


0.02 


0.29 


± 


0.01 


1.21 


+ 


0.05 


2.99 ±0.13 


0.50 


± 


0.07 


2.00 


+ 


0.00 


52.59 ± 0.38 


22 


F606W 


18.95 


± 


0.00 


0.77 


± 


0.02 


3.27 


± 


0.09 


3.78 ±0.13 


0.84 


± 


0.05 


3.00 


+ 


0.00 


46.77 ± 0.30 




F814W 


18.23 


± 


0.00 


0.83 


± 


0.03 


3.48 


± 


0.11 


3.94 ±0.16 


0.84 


± 


0.05 


3.00 


+ 


0.00 


54.90 ± 0.43 



Note. - Table[5]is published in its entirety in the electronic edition of the Journal. 
a Identification number in Paper I 
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